State and dynamical parameter estimation for open quantum systems 



O 
O 
(N 

Ph. 
< 



Jay Gambctta and H. M. Wiseman* 
School of Science, Griffith University, Brisbane 4^1^, Australia 

Following the evolution of an open quantum system requires full knowledge of its dynamics. In this 
paper we consider open quantum systems for which the Hamiltonian is "uncertain". In particular, 
we treat in detail a simple system similar to that considered by Mabuchi [Quant. Semiclass. Opt. 
8, 1103 (1996)]: a radiatively damped atom driven by an unknown Rabi frequency (as would 
occur for an atom at an unknown point in a standing light wave). By measuring the environment 
of the system, knowledge about the system state, and about the uncertain dynamical parameter, 
can be acquired. We find that these two sorts of knowledge acquisition (quantified by the posterior 
distribution for il, and the conditional purity of the system, respectively) are quite distinct processes, 
which are not strongly correlated. Also, the quality and quantity of knowledge gain depend strongly 
on the type of monitoring scheme. We compare five different detection schemes (direct, adaptive, 
homodyne of the x quadrature, homodyne of the y quadrature, and heterodyne) using four different 
measures of the knowledge gain (Shannon information about fi, variance in fi, long-time system 
purity, and short-time system purity). 

PACS numbers: 03.65.Yz, 42.50.Lc, 03.65.Wj 
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I. INTRODUCTION 

Quantum parameter estimation is a well-established 
area [1, 2], which is usually formulated as follows. A 
known quantum state enters an apparatus that performs 
an operation on the state. The operation, which is usu- 
ally unitary but need not be [3, 4], is parameterized by 
one or more unknown parameters. The goal is to esti- 
mate these parameters by making a measurement on the 
(unknown) output state. Except in special cases, it is 
not possible precisely to find out the unknown param- 
eters from a measurement on a single system. Rather, 
the operation and measurement must be performed re- 
peatedly, on a sequence of identically prepared quantum 
systems. 

There is a trivial sense in which it is possible to ob- 
tain complete information about the unknown parame- 
ters from a single system. That is by taking the out- 
put state after the measurement, and using it as the 
next input state, having perhaps transformed it first. If 
the transformation required is as difficult as preparing 
a new system from scratch, then there is nothing to be 
gained by reusing the same system. However, this sce- 
nario of repeated measurements on a single system is use- 
ful pcdagogically to make the transition to continuously 
monitored systems with unknown dynamical parameters. 
This transition is made by considering the limit where the 
unknown transformation is infinitesimally different from 
the identity, and the repeat time is infinitesimal. 

To the best of our knowledge, a theoretical treat- 
ment of estimating an unknown dynamical parameter 
by continuous observation of a system was first done by 
Mabuchi [5] . His system was a two- level atom coupled to 
a classically driven electromagnetic field mode in a cavity. 
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The unknown parameter was the position of the atom. 
This is a dynamical parameter because it determines the 
strength of the coupling between the atom and field (the 
Rabi frequency). The continuous monitoring considered 
was counting the photons that escape through one of 
the cavity mirrors. Mabuchi used Bayesian statistics to 
determine the posterior probability distribution for the 
Rabi frequency. This represents the knowledge the ex- 
perimenter would have about the Rabi frequency given a 
particular (typical) measurement record. The measure- 
ment is continuous in time (monitoring) because in any 
instant of time a photon may or may not be detected. 

In this paper we arc concerned with the same question, 
namely how would an experimenter gain knowledge of 
an unknown dynamical parameter from the measurement 
record resulting from monitoring the system. We even 
choose a similar (but even simpler) quantum system to 
that of Ref. [5], namely an atom driven by a classical 
field of unknown Rabi frequency. However, our analysis 
goes beyond, and has additional aims to, that of Ref. [5] 
(although we should note that extensions similar to the 
first three outlined below were suggested in a footnote of 
that work.) 

First, we consider the entire ensemble of possible mea- 
surement records and parameter values, rather than just 
one (typical) measurement record from one parameter 
value. 

Second, we quantitatively characterize this ensemble 
by calculating the average information gained (in bits) 
by the measurement, as a function of time. 

Third, we consider different ensembles resulting from 
different measurement schemes on the system. We em- 
phasize that the choice of measurement scheme does not 
affect the evolution of the system on average. That is, for 
all measurement schemes, averaging over the possible re- 
sults and the possible values of the Rabi frequency yields 
the same equation of motion for the system state. Physi- 
cally, this is because the average behaviour of the system 
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is determined by its immediate environment, whereas the 
different measurement schemes are effected by detecting 
the light emitted by the system in different ways. How- 
ever, the different measurement schemes give very dif- 
ferent typical posterior distributions, and very different 
rates of information gain. 

Fourth, and perhaps most distinctively, we consider 
not just the estimation of the unknown parameter, but 
also the estimation of the state of the system conditioned 
on the measurement results [6]. We do this using the 
same Bayesian method as for the parameter estimation. 
In this respect, our work could be seen as an extension of 
quantum trajectory theory [7] to systems with unknown 
dynamical parameters. Quantum trajectory theory is 
simply the application of quantum measurement theory 
to continuous monitoring of open quantum systems, most 
usually optical systems subject to photodetection [8] . 

If the dynamical parameters for an open quantum sys- 
tem are known then conditioning the system on efficient 
detection of its emissions is guaranteed to monotoni- 
cally increase its average purity in time, as information 
is gained about the system. But if dynamical param- 
eters are not known then the average purity may de- 
crease, as the different possible evolutions are summed in- 
coherently. On the other hand, the measurement record 
also contains information about these parameters, so that 
these parameters become better defined over time. Hence 
one might expect that the system will eventually become 
pure anyway. 

It is one of the main results of this paper that this 
expectation is not met. For our system there are some 
monitoring schemes for which the parameter never be- 
comes sufficiently well known for the system state to be- 
come pure. However, there is no simple correlation be- 
tween the information gained about the parameter (the 
Rabi frequency) and the final purity of the system (the 
atom). One monitoring scheme yields almost no param- 
eter information, yet produces, on average, a much purer 
final system state than do other schemes that yield large 
amounts of parameter information. Moreover, the rates 
at which the system state purifies is, for some monitoring 
schemes, tied to the rate of parameter information gain, 
while for other monitoring schemes it is much faster than 
that. These results can be understood only from an ap- 
preciation of the conditional dynamics induced by the 
different detection schemes. 

The remainder of this paper is organized as follows. 
In Sec. II we present the general formalism for state and 
dynamical parameter estimation by monitoring a single 
system. We also explain how the parameter information 
gained is quantified. In Sec. HI we introduce the sys- 
tem to which we apply our formalism, a two-level atom, 
driven by an unknown Rabi frequency, and monitored by 
having its fluorescence detected. Sec. IV contains the re- 
sults of our numerical simulations of the relevant ensem- 
ble averages for five different detection schemes: direct, 
the adaptive scheme of Wiseman and Toombcs [9] , homo- 
dyne of the X quadrature, homodyne of the y quadrature. 



and heterodyne. Sec. V concludes. 



II. GENERAL FORMALISM 
A. Quantum trajectories 

It is well known that quantum trajectories can be used 
to describe the evolution of a continuously monitored 
open system [8] . Since here we are continuously monitor- 
ing an open system with an unknown dynamical param- 
eter, we begin by giving a brief outline of the standard 
quantum trajectory theory. 

A good place to start is with the measurement formal- 
ism for open systems [10, 11]. An open system is simply 
a quantum system that interacts with its environment 
(usually called a bath). This interaction, like all quan- 
tum interactions, generally entangles the system and the 
bath. If we initially have states |V'(^o)) and |m(to)) for 
the system and bath respectively, and let these states en- 
tangle by [/(to + r), a unitary operator which includes 
both the bath-system coupling and the system dynam- 
ics. An instantaneous rank-one projective measurement 
on the bath will result in the state after the measurement 
being 

^^^,u^T^\ k) {r\U(to + T)\mit„)) mo)) 

\r) \iPr[to + T}} = -== , 

(2.1) 

where P(r) is the probability of getting the result r. 
Eq. (2.1) shows that after the measurement the system 
and the bath are disentangled, so it is not necessary to 
continue to describe the bath in our treatment of the 
measurement. This allows Eq. (2.1) to be reduced to 

I , f, ^^^^ MrjT) |V>(t)) 

\A{to + T)} = == , (2.2) 

where Mr{T) = (r| U{ta + T) |m(to)) is called the mea- 
surement operator and has the feature of collapsing the 
observer's knowledge of the system into a state that is 
consistent with the result r. Mr{T) is still an operator 
for the system as U (to + T) is an operator on the ten- 
sor product Hilbert space for system and the bath. It is 
important to note that this measurement operator is not 
necessarily a projector in the system Hilbert space. 
The probability P(r) is given by 

P(r)=Tr[F,(r)|V(to))(V'(io)|], (2.3) 
where Fr is called the effect and is defined as 

Fr{T) = Ml{T)Mr{T). (2.4) 
The complete set of effects must sum to one: 

Y,Fr{t)^l. (2.5) 
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The above formalism for measurement only considers 
pure states, but to take into account initially mixed states 
Eq. (2.2) can be rewritten in terms of the state matrix. 
The state after the measurement is then 



Pr{ta + T) = Mr{T)pm)Ml(T)/V{r). 



(2.6) 



Here Eq. (2.6) describes the state conditioned on the re- 
sult r and is referred to as an unraveling of the average 
post-measurement state p(to + T). That is, the weighted 
mean of all the possible conditioned states for one unrav- 
eling is equal to the average state: 

p(io + T) = E[p,(to + T)]=Y^ P(r)p.(to + T). (2.7) 



It should be noted that an average state has more then 
one unraveling. The different unravelings correspond to 
different sets of measurement operators, arising from dif- 
ferent sets of environment projectors \r) (r| in Eq. (2.1). 

As mentioned earlier, quantum trajectories arises when 
this measurement formalism is applied to a continuously 
monitored open system [8]. In continuous monitoring, 
repeated measurements of duration 5t arc performed on 
the state. This results in the state being conditioned on 
a record I[o.t), which is a string containing the results rfe 
of each measurement. Here the subscript k refers to a 
measurement at time tk = k5t, with ~ Q. Using this 
I[o (), the conditioned state at time t can be written as 



Pi(t) -pi(t)/P(I[o,t)), 



(2.8) 



where pi{t) is an unnormalized state conditioned on I[o,t) 
and is equal to 

h{t) 

= M,,M,,_, . . . Mr.pmil ■ ■ ■ Ml_Ml. (2.9) 

The probability of obtaining this record is 

P(I[o,t)) = P(rfc)P(rfc_i) . . . P(ri) - Tr[pi(t)]. (2.10) 

To completely achieve continuous monitoring we let 
the time step between measurements, 5t, tend towards 
the infinitesimal interval dt. In doing this, Eq. (2.8) de- 
fines a stochastic master equation (SME) , with its ensem- 
ble average reproducing the usual deterministic master 
equation. That is, 

p{t) = ^p(Vo)pi(t) = 5^Pi(0 

l[0,t) l[0,t) 

= ,^^0 ^ Mr,„,...Mr,p{Q)Ml---Ml^^^ 
= lim (1 + C5tf'^* p{Q) = exp(£t)p(0), (2.11) 

5t — *0 

where for arbitrary p, C is the Liouvillian superoperator 
defined as Cp = \\mst-*o{Yl,r M^pMl - p)/5t. 



B. Quantum trajectories with an unknown 
parameter 

We now consider the situation where there is an un- 
known dynamical parameter A in £, and hence in the 
measurement operators Mr- This is done by simply not- 
ing that for each A there will be a conditioned state. This 
gives a doubly conditioned state of the form 



Pi,a(<) =Pi,A(i)/P(I[o,ol^). 



(2.12) 



where P(I[o_t)|A) is the probability of getting I[o,t) given 
A. It is obtained by 



P(I[o,*)|A) = Tr[pi,;,(i)]. 



(2.13) 



We wish to determine the posterior probability distri- 
bution P(A|I[o^t)) of A, given I[o.t)- This can be achieved 
using a Bayesian inference formula [12]. 



P(A|I[o,t)) 



P(I[o,,)|A)Po(A) 
/P(I[o.t)|A)Po(A)dA' 



(2.14) 



where Po(A) is the prior distribution for A. For a "good 
measurement" of A, as time increases, we would expect 
this prior distribution to converge to a (5-distribution. 

Theoretically, Eq. (2.14) is complete for determining 
P(A|I[o,t))- However, in general P(I[o^t)|A) is very small 
and in numerical simulations it will incur large computer 
roundoff errors. The small magnitude of P(I[o,t) | A) is due 
to the many possible trajectories the system could follow. 

To overcome this problem, linear quantum trajectories 
[13] were used. Linear quantum trajectories arise if we 
assume an ostensible distribution for the result r, A(r) 
[8]. These A(r) arc independent of A and the only con- 
dition they must satisfy is that they add to one. With 
these ostensible probabilities, the linear stochastic mas- 
ter equation (LSME) is derived from [8] 

Pi,A(i) -/5i,A(t)/A(I[o,t)), (2.15) 

where the ostensible probability for getting I[o,t) is 

A(I[o,f)) - A(rfe)A(rfe_i) . . . A(ri). (2.16) 

The actual probabihty of getting I[o^i) is [8] 

P(I[o.t)|A)=A(%,))Tr[pi,A(t)]. (2.17) 

Substituting P(I[o.t)|A) into Eq. (2.14) we obtain 

Tr[pi,A(t)]Po(A) 



P(A|I[o,t)) - 



/ Tr[pi,A(<)]Po(A)dA" 



(2.18) 



From Eq. (2.18) wc see that to calculate P(I[o_t)|A); 
the norm of the linear conditioned state [Eq. (2.15)] is 
needed. The order of magnitude of this norm is de- 
pendent on the ostensible probability we chose. By 
Eq. (2.17), if A(I[Q ()) is chosen to be of the same order 
as the true probability, this norm will be of order unity. 
This avoids the problem of large computer roundoff error. 
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C. Quantifying the information gained 

One of the main aims of this paper is to classify the 
information gained about the unknown parameter. The 
posterior probabihty calculated by Eq. (2.18) contains all 
the information about A for a particular record. However 
the question remains, how can this information be quan- 
tified? Two measures were investigated. The first is the 
variance: 



Vi = I P{X\I[o,t))X'dX-(^l P(A|I[o, 



t^)Xd\j . (2.19) 
The second is the information gain, A/j defined as [14] 
A/i = J P(A|I[o,t))log2P(A|I[o,t))dA 

Po(A)log2Po(A)dA. (2.20) 



This measures the number of bits of information gained 
by the observer about the parameter A. It can be thought 
of as the negative change in entropy of A. The great- 
est information gain corresponds to the transition from 
a flat (most disordered) distribution to a peaked (most 
ordered) distribution. 

These parameters give an indication of the quality of 
knowledge gained by an observer, for a particular run of 
the experiment. To characterize a particular measure- 
ment scheme, it is necessary to calculate the ensemble 
averages of Vi and A/j, which we denote as V and A/. 
The ensemble average of a parameter Ai is defined as 

A =E[Ai] = ^AiP(I[o.t)) 

1(0. t) 



J AiP(I[o,t)|A)Po(A)dA. (2.21) 



I[o,t 

Numerically, this is done by picking a true A, Atruci 
randomly from Po(A), and then simulating a quantum 
trajectory for this Atruo, yielding I[o,t)- This gives a typ- 
ical record as would be obtained experimentally. This 
I[o,t) is then used to calculate Tr[pi^\{t)] for all A's in the 
range of Pq. This allows the calculation of P(A|I[o_t)), 
with this probability the parameter of interest Aj can be 
calculated. By storing this value and repeating the above 
procedure n » 1 times, the ensemble average A of Ai is 
obtained. 



D. Best estimate of conditioned state 

Another aim of this paper was to determine the best 
estimate of the state given the knowledge we have ob- 
tained from a measurement. In Eq. (2.12) we defined the 
doubly conditioned state that arose when the state was 
conditioned on both 1[o^t) a-nd A. From Eq. (2.12) there 
are two best estimate states that can be calculated. They 



are p\ and pi and can be interpreted as the best estimate 
state, when A or I[o,t) is known respectively. They are de- 
fined as follows 



/ pi.A(t)Po(A)dA 
/P(I[o,i)|A)Po(A)dA- 



(2.22) 
(2.23) 



It should be noted that the average of each of these states 
will give the same average state p{t). 

Equation (2.22) describes the best estimate state that 
arises when the dynamical parameter is known and the 
record is not (i.e. a non-monitored system). This obeys 
master equation px = C\p\. Of more interest to us is the 
best estimate state described by Eq. (2.23), which is the 
state conditioned on some observed record I[o,t), when 
the true value of A is unknown. 

In calculating pi, if we use Eq. (2.23), we again run 
into the problem that the magnitude of pi,\{t -f dt) will 
typically be very small. Again this is overcome by using 
linear quantum trajectories, replacing Eq. (2.23) by 



Pi{t) 



J pi,A(^)Po(A)rfA 
/Tr[pi^A(i)]Po(A)dA- 



(2.24) 



To quantify the information gained about the state, 
the purity (pi) can be determined. 



PI = Tr[pi(<)2 



(2.25) 



The ensemble average purity {p ~ E[pi]) will give us an 
indication of how well the measurement scheme is at pro- 
ducing pure states. One might expect that a highp would 
correspond to a high AI. However it will be seen that 
this is not true. 



III. THE SYSTEM 

The system we are considering is a classically driven 
two level atom, immersed in the vacuum. With no mon- 
itoring of the vacuum field, the average state evolution 
when all the dynamical parameters are known is given 
by the master equation. The Lindblad form [15] of the 
master equation for the TLA, in the interaction picture 
(with respect to the free evolution of the atom) is [16] 



^V[a]p{t) ^ Cnp{t). (3.1) 



Here is the Rabi frequency, 7 is the spontaneous emis- 
sion rate, a is the lowering operator, is the usual Pauli 
matrix and V is the superoperator that represents damp- 
ing of the system into the environment. It is defined as 
[17] 



(3.2) 
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The solution of this equation can be described by the 
Bloch vectors {x,y, z), with p written as 



^(1 



yay 



The purity p is equal to 



(3.3) 



(3.4) 



Using this Bloch representation the solution of 
Eq. (3.1) is a state that rotates about the x-axis at 
frequency Q, with damping in all variables towards the 
steady state value of 



0, 



Vs: 



-7 



21^2 + 



2172 + ^2 



(3.5) 



The most obvious choice for the unknown dynamical 
parameter is ft, as indicated by the subscript in Cn in 
Eq. (3.1). This can be physically motivated as follows: 
if we placed a laser-cooled atom (with no center-of-mass 
motion) in a classical standing field, then the il it would 
experience is 



(3.6) 



where k is the wavevector for the classical field and x is 
the position of center of mass of the atom. We assume 
that the placement of the atom in the field is not biased 
in any way. That is, in one wavelength (A) of the field 
the atom position distribution is given by Po(a;) = 1/A. 
Using Eq. (3.6), Po(a^) can be transformed into a proba- 
bility distribution in CI space. 



Pom 



1 



2 - (12 
max ^ ' 



(3.7) 



This is the prior distribution for il, that will be used in 
the rest of this paper, with fimax = IO7. Along with this 
prior distribution the initial condition that we will use 
for our simulations, unless otherwise stated, is pi,n(0) 
satisfying 



CnPi.niO) - 0. 



(3.8) 



That is, we will assume the initial state is the steady 
state of the general master equation Eq. (3.1). 



IV. RESULTS 

The results of this paper are broken down into five 
subsections, each corresponding to one of the five mea- 
surement schemes investigated. 

A. Direct Detection 

The first measurement scheme investigated was direct 
detection. This involves the detection of all the fluores- 
cence emitted by the atom as shown in Fig. 1. Con- 
tinuous monitoring with this detection scheme will yield 



Photodetector 



FIG. 1: A schematic for direct detection. The atom is placed 
at the focus of a parabolic mirror so that all the fluorescence 
emitted by the atom is detected by the photodetector. 



either one of two results for each interval dt, a detection 
(labelled by a 1) or no detection (labelled by a 0). Thus 



will be a string of O's and I's. The measurement 



operators for each of these results are [8] 
Mi{dt) = \/dff<7, 
Moidt) = 1 - (i^a, + ^a^a ] dt. 



(4.1) 
(4.2) 



It can be shown that these measurement operators sat- 
isfy the completeness condition, Eq. (2.5). Using these 
measurement operators and Eq. (2.8), a SME for direct 
detection can written as 



dp\.n = dN [t)Q{^ dt-i a\pY.n - dtn[i—crx + ^crVJpLo, 

(4.3) 

where G and Ti. are the nonlinear superoperators defined 
for arbitrary a and p by 



g[a]p 
n[a]p 



apa' 



Tr [apa^"] 



P, 



(4.4) 



ap + pa' — Tr[ap -|- pa']p- (4.5) 



In Eq. (4.3), the variable dN is a stochastic increment 
that equals one if there is a detection in the interval dt 
and equals otherwise. Formally, dN is defined by 



dN{tf = dN{t), 
E[dN{t)] = P{l) = dt-f{a^(j). 



(4.6) 
(4.7) 



By averaging Eq. (4.3) and using Eq. (4.7), it is easily 
seen that the SME is an unraveling of the general master 
equation, Eq. (3.1). A typical trajectory of this SME is 
shown in Fig. 2 (solid line), for fl = 5-f. It is observed 
that the x component is zero, and the y and z oscillate 
in quadrature. This can be understood physically as the 
state is dominated by the ^a^/^ Hamiltonian, with de- 
tections occurring stochastically according to Eq. (4.7). 
After each detection the state collapses to the ground 
state {x = 0, y = 0, and z = —1). 

To consider the case when Ct is unknown, a LSME had 
to be developed. Using the direct detection measurement 
operators and Eq. (2.15), with A(r) defined as 



A(l) = edt = 1 - A(0), 



(4.8) 
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FIG. 2: The best estimate states when Q is known (Solid 
hne) and unknown (dotted hne) for iltmo = 57 when direct 
detection is used. Time is measured in units of 7~^, x, y, 
and z are the Bloch vector components and the purity p — 
^{l + x^ +y'^ + z^). The initial states are the steady state for 
the known case and the average steady state for the unknown 
case. 

where e is an arbitrary parameter, the LSME is 

^dtn[i^<7, + |(tV - (4.9) 
The g and Ti. Hnear superoperators are defined as 

G[a]p = ^-P- (4-10) 

n[a]p = ap + pa''. (4.11) 

To obtain the general master equation from Eq. (4.9), 
E[diV] = A(l) has to be used. However, to determine the 
parameters of interest to us, namely pi{t) and P(r2|I[o,t)), 
Eq. (4.9) is numerically simulated for all possible fl in 
Po(ri) with dN specified by I[o,t)- This record would 
ideally be obtained experimentally but for the purpose 
of this paper it is calculated by numerically evaluating 
Eq. (4.3) for a known fl, which we will refer to as Sltruo- 
This I[o,t) is then used in Eq. (4.9) to generate Tr[pi,a] 
for all the f2's between — f2max and f2max- Then with 
Eq. (2.24) and Eq. (2.18) one can obtain both pi{t) and 

For a I[o,t) based on fltmc ~ 57 the best estimate state 
and the posterior distribution where calculated and are 
shown in Fig. 2 and 3 respectively. It is observed that, 
in contrast to the known case, the best estimate of y is 
identically zero. This is because positive and negative 
are initially equally likely, so that i/ss in Eq. (3.5) averages 
to zero. Moreover, the sign of fJ is not determinable by 
this measurement scheme, because the rate of detections 
depends only on z, which is independent of the sign of Q. 

Another difference apparent with the unknown CI case 
is that z oscillates with a different frequency to the known 




FIG. 3: A plot of a typical P(f2|I[o,t)) when Sltruo = 5 for the 
direct detection scheme, fl is measured in units of 7 and time 
is measured in units of 7~^. 



case, in this case a faster frequency [since Po(f^) is 
peaked at the end points = S^max = IO7]. However as 
time increases its frequency tends to that of the known 
case. This is due to the fact that for direct detection 
the rate of detections is dependent on the magnitude of 
f2, so as time goes on one would expect to gain more 
information about the magnitude of Q. 

These interpretations of the conditioned dynamics are 
confirmed in Fig. 3. With increasing time, the posterior 
distribution localizes at ztfitrue. The mean is always zero 
and thus is not an unbiased estimator of fl. The reason 
that the magnitude is determinable and the sign is not, 
can be formulated as follows. In the Bloch representation 
of Eq. (4.9), with the transformation y —y, fl —fl 
the equations stay invariant. Since this transformation 
changes the direction of rotation around the x-axis, we 
will call it the rotation transformation. 

With an indeterminable direction of rotation and this 
measurement scheme, it can be seen that the best esti- 
mate state will never become more pure than a state that 
is a mixture of two states that rotate is opposite direc- 
tions around the a; = great circle of the Bloch sphere. 
Thus the best estimate state oscillates up and down the 
z-axis of the Bloch sphere. 

We turn now to quantifying the measurement scheme's 
ability to gain knowledge, by numerically determining the 
ensemble average purity, V and A/. These ensemble av- 
erages were calculated for ^Itruc's weighted on the prior 
distribution, Eq. (3.7). These numerical simulations are 
depicted in Fig. 4 for two initial states; one is the steady 
state (solid line) and the other is the ground state (dot- 
ted line). It is observed that in both cases the average 
purity of the state never attains one, with the purity in 
the second case initially decreasing from one. The long 
time purity (~ 0.75) is due the best estimate being a mix- 
ture of two states as explained above. This figure can be 
obtained analytically, if we make the follow two assump- 
tions. The first is that Cltmc 3> 7. This is valid as Po{^) 
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FIG. 4: The ensemble average (n = 1000) of the purity, vari- 
ance and A7 when direct detection is used, for two initial 
states, the steady state (solid) and ground state (dotted). 
Time is measured in units of 7^^. 



from which fitruc is drawn is peaked at ztrimax; and in 
our calculations llmax 2> 7. The second assumption is 
that in the long time limit the posterior distribution lo- 
calizes on zbOtrue, which is what is seen in Fig. 3. With 
these two assumption the long-time best estimate state 
in Bloch representation will be 



0, 



y = 0, 



'C0Sr2truc(i - ilast), (4.12) 



where tiast is the time of the last jump, which is typically 
more than one Rabi cycle before t. With this state the 
average purity (for the long time limit) can be estimated 
as 



2TT/n 



1 



2 4 



(4.13) 



From Fig. 4, it is also observed that the simulated en- 
semble average variance V is approximately constant for 
all time. In fact, given that the no information about 
the sign of O is determinable, and that the initial distri- 
bution Po(fi) is symmetric, it is easy to prove that V is 
exactly constant. 

For the third parameter A/, it is observed that, on av- 
erage, direct detection yields information about as time 
increases, for both initial states. It is observed that the 
initial slope of A/ is zero for the ground state, while it is 
non-zero for the initial steady state case. The initial flat- 
ness in the first case is due to the fact that if the system 
starts in the ground, the rate of detections (proportional 
to the excited state component) scales as (iltruci)^, and 
with out any detections it would not be possible to gain 
any information. By contrast, for the steady state case 
there will be some excited state fraction (depending on 
il) and thus a finite detection rate even at t = 0. Fig. 4 
also show that, after the initial flatness, the AJ in the 
first case rapidly overtakes that in the second case. This 
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FIG. 5: A schematic for adaptive detection. The fluorescence 
emitted by the atom is coherently mixed with a weak local 
oscillator (LO) via a low reflectivity beam splitter (LRBS). 
The electro-optic modulator (EOM) reverses the amplitude 
of the LO every time the photodetector fires. 



jump in AI occurs at roughly t — l/r^max; which is when 
one would expect a significant excited state fraction to 
have developed (Recall that Po(f^) is sharply peaked at 



— ztr^rnax)- 



B. Adaptive Detection 



The second measurement scheme investigated was the 
adaptive scheme of Wiseman and Toombes [9]. For a 
known Jl, this measurement scheme is designed to keep 
the atom jumping between two fixed states. For large, 
these fixed states turn out to be close to ax eigenstates. 
This two-state jumping is achieved by coherently mix- 
ing the fluorescence emitted from the atom with a weak 
local oscillator (LO) via a low-reflectance beam splitter 
(see Fig. 5). The reflected amplitude ^ of the local oscil- 
lator is switched between zbi^/7 each time a detection is 
registered by the photodetector. 

For this detection scheme the measurement operators 
arc [9] 



(4.14) 



Mi{dt) 

Moidt) 

(4.15) 

These measurement operators result in a SME of the form 
dpi,n = dN{t)g[y/dt^ {a + - dm[C]pi,n, (4.16) 

where 

C^i^'yx + ^<y^cT + fija+^. (4.17) 

Using the same ostensible distribution A(r) as in direct 
detection, the LSME is 

dpi.n = dNit)g[^dfi (a + p)]pi,n ~ dtn[C]pi,n, (4.18) 
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FIG. 6: The best estimate states when adaptive detection is 
used. Details are as in Fig. 2. 



where 



c 



i — a.r 
2 



^ t 
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(4.19) 



Figure 6 shows the best estimate state for a known 
(sohd) and unknown U, (dotted), with ritiuc = 57. It 
is observed that with the known case after the initial 
transients, the state jumps between the two fixed states 

[9] 



202 _^ ^2 



y = 



2O7 
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202 + ^2 ' 



202 



■7" 
(4.20) 

For the unknown case the y component averages to 
zero, and the x and z components both appear to be 
sHghtly different to the known case. 

Similarly to the direct detection case, a better under- 
standing of this state can be obtained by considering 
P(r2|I[o,t)). This is shown in Fig. 7 and it can be seen 
that as time increases under this adaptive measurement, 
the typical posterior probability distribution P(f2|I[o^t)) 
scarcely changes from Po(ri). This is not unexpected, as 
for this detection scheme it can be shown that at steady 
state the jumps are Poissonian, with rate 7/4. That is, 
the jumps are independent of [9] and hence yield no in- 
formation about it. Since P(r2|I[o.t)) ~ Po(r2), we can use 
this approximation to obtain analytically an indication 
of the best estimate state by solving Eq. (2.23). For this 
detection scheme this is simply the mean of Eq. (4.20) 
under the distribution Po(r2). This gives 



x = T(l- 



7 



y = 0, z = 



-7 



2 

max 



+ 72 

(4.21) 

Comparing this with the numerical simulation it is ob- 
served that they agree very well. 

To quantify this detection scheme, the ensemble av- 
erage of the variance, purity and AJ were numerically 
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FIG. 7: A plot of P(f2|I[o.t)) for the adaptive scheme. Details 
are as in Fig. 3. 




FIG. 8: The ensemble average (n = 1000) of the purity, vari- 
ance and A/ when the adaptive detection technique was used. 
Note for A/ the scale has been change when compared to 
Fig. 4. 



calculated and are shown in Fig. 8. The purity rapidly 
becomes, and remains, relatively high. This is because 
the best estimate state of Eq. (4.21) is the same no mat- 
ter what f2truc is chosen. For fimax = 107 the numeri- 
cal value of the stationary purity is 0.934 and by using 
Eq. (4.21) an analytical value of the purity can be ob- 
tained, 



7 



7 



r 



2 a, 



(4.22) 



For ilinax = IO7 this gives a value of 0.934, which is equal 
to the numerical value. 

Since this state has a high purity one might expect 
that the unknown parameter must also be well defined. 
However this is not true as already discussed. This lack 
of knowledge about Q. is seen in Fig. 8. Like direct detec- 
tion, the sign of cannot be determined so the average 
variance remains precisely constant. However unlike di- 
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rect detection, the information gain is bounded, with a 
maximum AI of less than 0.06 bits. 

An interesting point to note is this scheme would be 
well suited to estimating 7 (if there were some uncer- 
tainty in that parameter) even if Q was also uncertain. 
That is because the detection rate is proportional to 7, 
almost independent of Q. Of course this would require 
the local oscillator amplitude to be adjusted from an ini- 
tial guess according to the best estimate of 7. 

C. Homodyne x Detection 

To perform a homodyne detection experiment, a simi- 
lar arrangement to the adaptive scheme is used. That is, 
the output flux from the atom is mixed with a resonantly 
tuned LO by a beam splitter (see Fig. 9). However in this 
scheme there is no feedback and the amplitude /? of the 
LO is assumed to be infinite {f3 00). Because of this, 
there will be many detections in the interval dt. Each de- 
tection causes only an infinitesimal change in the system 
state, so the evolution of the system can be described by 
a diffusive SME. In each dt there will be a continuous 
current I registered in I[o,t) rather than a detection or 
no-detection. Since / is a continuous variable we can de- 
fine a measurement operator, Mj, a continuous function 
of /, to represent this measurement scheme. 

Mi = v/t7[1 - {i^a, + ^a^a - ^ ae"''' I)dt]. (4.23) 

Here $ is the phase of the local oscillator and 

Tidl = —p^=e-^^^'^^dl, (4.24) 

is a Gaussian probability measure. It is easily shown that 
this continuous measurement operator satisfies the com- 
pleteness condition, Eq. (2.5), where the sum is replaced 
by an integral over / between ±00. 

With this continuous measurement operator the SME 
in the Ito form is [18] 

dpi.n = Cnpi,ndt + y/jH[ae~^'^]pi 

X {Idt - V7 Tr[(Te-'*pi,n + pi,aate**]di), 

(4.25) 

where / is the current element for the interval dt and 
is equal to the difference between the number of detec- 
tions at the two photodiodcs divide by the intensity of 
the field. By using Eq. (4.23) and Eq. (2.3) the proba- 
bility of getting / for the interval dt can be calculated. 
This gives a Gaussian distribution with a mean equal to 
(cre~'* -I- cr^e'*) and a variance of dt^^. Thus, / will 
be a Gaussian random variable (GRV) of the form 

/ ^ V7Tr[ae-'*pi + pia^e**] + ^(t), (4.26) 

where £^{t) = dW{t)/dt represents Gaussian white noise, 
and is formally defined as [21] 

E[^(i)] = o, mmn^sit-t'). (4.2?) 




FIG. 9: A schematic for the three detection schemes, homo- 
dyne of the X quadrature, homodyne of the y quadrature and 
heterodyne. For the homodyne schemes the LO is resonantly 
tuned to the atomic frequency, with a phase of zero and n/2 
for the X and y schemes respectively, whereas for the hetero- 
dyne it is detuned by an amount A. 

For the LSME we take the ostensible probability for 
the current to be equal to that which would arise from 
the LO alone. This results in A(/) = T/ so that I is 
ostensibly a GRV with mean zero and variance dt~^, like 
C(t). The LSME in Ito form is [18] 

dpi,n = Cnpi.ndt + ^n[ae-"^]pynldt. (4.28) 

It can be seen that both the LSME and the SME re- 
duce to Eq. (3.1) when the ensemble average is taken. 
Similarly to the previous schemes, to determine an un- 
known n, I[o_t) is generated by the SME for a preset fl, 
^true, which may then be "forgotten". The LSME is 
then used to generate both pi{t) and P(r2|I[o,t)) for the 
predetermined record I[o,t)- 

For homodyne x quadrature measurement, the $ of the 
LO is set to zero (as a; = (cr + (t^)). With this phase and 
f^true = 57, the best estimate state for a known and un- 
known Q are shown in Fig. 10. It is observed that for the 
known case, the state seems to localize itself relatively 
fast into pure states that have a large x contribution, and 
small oscillations in the y and z directions. By contrast, 
when is unknown, the best estimate state still contains 
a large x contribution, but the y is strictly zero and the 
amplitude of the z oscillations is reduced. As in the pre- 
vious cases, this zero y component can be understood 
by considering P(ri|I[Q j)), shown in Fig. 11. It is seen 
that, like direct detection, this measurement scheme has 
an even posterior distribution that localizes at ±r2truo- 
This is again due to the stochastic Bloch equations be- 
ing invariant under the previously considered rotational 
transformation. However, the rate at which this localiza- 
tion occurs is much slower than under direct detection. 

The slower rate of information gain is confirmed with 
the calculation of the ensemble average of A/, shown in 
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FIG. 10: The best estimate states for the homodyne x scheme. 
Details are as in Fig. 2. 



FIG. 12: The ensemble average (n = 500) of the purity, vari- 
ance and A/ when homodyne x was used. Time is measured 
in units of 7~^. 
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FIG. 11: A plot of P(f2|I[o,t)) for the homodyne x scheme. 
Details are as in Fig. 3. 



Fig. 12. It is seen that within 507^^ units of time, A/ for 
homodyne x is about half that of direct detection. Phys- 
ically this comes about because, for the system we are 
investigating, the underlying dynamics cause the states 
to rotate around the x-axis with frequency fitrue- The 
measurement scheme tends to produce states oriented 
mainly in the ±a; directions. This can be understood 
from the measurement effect F/, which, using Eq. (2.4), 
can be shown to be 



F,dl 



1 



^27r/dt 



le 2 



dl. 



(4.29) 



This effect is a Gaussian with a mean equal to the Ux 
quadrature operator and variance dt-^. Thus, it IS an 
unsharp measurement of ax- Thus, for a measurement 
scheme that makes the conditioned state mainly oriented 
in the ±x directions, one would expect that this state 
would be less affected by an unknown O than a state on 
the X = plane as produced by direct detection. Thus 
less information about 17 comes out of the measurement 



record. In Fig. 11 it is observed that the ensemble aver- 
age of the purity of this state increases quickly to about 
0.75, then increases only slowly afterwards. This quick 
increase is also a result of the state becoming predomi- 
nately ±a; oriented, (similarly to the adaptive detection 
scheme) and the slow increase is due to the slow increase 
in the knowledge of VL (similarly to direct detection) . As 
with direct detection, the system state will never become 
fully pure. This is due to the double peaks in P(r2|I[o,t)), 
which insures the y component of the state always aver- 
ages to zero. 



D. Homodyne y Detection 

Setting the $ of the local oscillator to 7r/2 allows mea- 
surement of the y quadrature (as y = (— zcr -I- icr^))- The 
best estimate states for the known and unknown f2 are 
shown in Fig. 13, for f2true = 57. It is seen that when O is 
known (solid) this measurement scheme makes the state 
coarsely rotate around the Bloch sphere with a purity of 
one. When Q, is unknown (dotted line), Fig. 13 shows 
that, unlike the previous schemes, the y component does 
not average to zero. As time increases the oscillations in 
the y and z components for the unknown SI case gradually 
converge to those for the known O case. This suggests 
that this scheme can determine ritruc- This is confirmed 
by the calculation of P(f7|I[o^t)) shown in Fig. 14. The 
ability of this scheme to distinguish the sign of 17 can be 
physically understood by considering the Bloch represen- 
tation of Eq. (4.28). These stochastic equations are not 
invariant under the rotation transformation. 

To understand how this scheme reduces the uncer- 
tainty in Jl, consider the effect for this measurement 
scheme 



Fidl = 



1 



yfhi/dt 



(4.30) 
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FIG. 13: The best estimate states for homodyne y measure- 
ment. Details are as in Fig. 2. 
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FIG. 14: A plot of P(n|I[o.t)) for homodyne y measurement. 
Details are as in Fig. 3. 

That is, Fi is an unsharp measurement of y. Now y is 
a variable that is directly affected by f2truo, and indeed 
the sign of y reverses if the sign of O reverses. Even 
though in each interval dt, y is measured unsharply, over 
time this detection scheme will result in a narrowing of 
our knowledge of Q, until infinite time where it would be 
fully known. This is further confirmed by the calculation 
of the ensemble averages of the three parameters, purity, 
V and AI (Fig. 15). It is observed that the purity of this 
state increases up to one, the V inD. reduces substantially 
in the 507~^ units of time and AI increases to a value 
larger than that for all other schemes. 

E. Heterodyne Detection 

The last detection scheme considered uses the hetero- 
dyne technique. This detection scheme uses the same 
arrangement as the homodyne (see Fig. 9), with the only 
difference being that the LO is now detuned from the 



FIG. 15: The ensemble average {n = 500) of the purity, vari- 
ance and AI when homodyne detection of the y quadrature 
was used. 

atom by an amount A. This effectively results in the LO 
having a time varying phase of At with respect to the 
driving field. Since the field amplitude is still assumed 
to be infinite as in the homodyne case, I[o,t) will com- 
prises of a string of real numbers /. However, by coarse- 
graining to obtain the Fourier components at w = ±A, a 
complex photocurrent is obtained [17]. The continuous 
set of measurement operators after the coarse graining 
approximation [Adt ^ 1 but jdt 1) are 

Mi = yT7[l - + l<j^a - ^ar)dt], (4.31) 

where 

Tjd^I = '^e-^'^'^'dH. (4.32) 

IT 

It is easily shown that these measurements operators sat- 
isfy the completeness condition, Eq. (2.5). To do this, one 
must integrate over the plane of the complex currents /. 

As with homodyne, the sample path for / can be ob- 
tained from calculating the probability of getting / in the 
interval dt. Doing this, one obtains 

I^[Vl{^)+m], (4.33) 

where C,{t) is a complex Gaussian white noise term, which 
is formally defined as [21] 

nman = nm] = o, (4.34) 

E[C(t)C(i')] = 5{t-t'). (4.35) 

Using the above measurement operators and 
Eq. (4.33), the heterodyne SME in Ito form is [18] 

dpisi. = \/l{(^Pi,n - (o-) pisi.){I*dt - (o"^) dt) 
+V7(pi,ofT^ - (a^) piM){Idt - V7 {(^) dt) 
+CnPiMdt. (4.36) 
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FIG. 16: The best estimate states, when heterodyne is used. 
Details are as in Fig. 2. 

For the LSME wc again assume that the ostensible 
probability is that due just to the LO, which results in 
a heterodyne current / with the same statistics as C(i)- 
With this complex current the ostensible probabihty A(/) 
is equal to T/. This gives a LSME in Ito form of [18] 

dpi,n = CnPi.ndt + y^ api^nl*dt + y^ pi,n(7^ Idt. (4.37) 

Using r^truc = 57, the best estimate state for known 
and unknown fl are shown in Fig. 16. It is observed that 
for a known fi, the state contains attributes of both the 
homodyne x and y measurement schemes. By this we 
mean that the state tends to have a distinct x compo- 
nents, whilst keeping the coarse rotations of the homo- 
dyne y scheme. This is not unexpected as heterodyne is 
equivalent to simultaneous homodyne x and y measure- 
ments, each of 50% efficiency [22]. In the unknown 
case it is observed that the y component does not aver- 
age to zero, suggesting that P(il|I[o,t)) focalizes to r^truc, 
which is confirmed by Fig. 17. However, the rate at which 
P(il|I[o,t)) converges to 6{^l — fltmo) is much slower than 
that of the homodyne y measurement. This is also il- 
lustrated in Fig. 18 as the ensemble average A/ is not 
as high. Fig. 18 also shows the ensemble average of the 
purity and from this figure it is seen that it contains sim- 
ilar properties of both the homodyne x and y schemes. 
In particular, it has an initial sharp increase, which is 
due the state obtaining a large x component (similar to 
the homodyne x scheme) and as time goes on the pu- 
rity increases to one due to the localization of P(ri|I[o_t)) 
(similar to homodyne y). 

V. DISCUSSION 

The results of this paper demonstrate that quantum 
parameter and state estimation for a continuously mon- 
itored open system is greatly affected by the measuring 
scheme. It was observed that as the measurement time 
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FIG. 17: A plot of P(f2|I[o.t)) for heterodyne detection. De- 
tails are the same as Fig. 3. 
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FIG. 18: The ensemble average (n = 250) of the purity, vari- 
ance and A7 for heterodyne detection. Time is measured in 
units of 7^""^. 

increased, some detection schemes had the ability of both 
reducing our uncertainty in the unknown dynamical pa- 
rameter, and producing a conditioned state of high pu- 
rity, whereas other schemes could only do one of these, or 
none (depending on how the uncertainty in the unknown 
parameter is quantified). We re-emphasize that all of 
the measurement schemes arise from the same coupling 
of the system to the environment; all that is different is 
how the environment is measured. 

The system we considered was a two-level atom with 
Hamiltonian flax/2, with spontaneous decay rate 7. The 
unknown dynamical parameter is Q, the Rabi frequency. 
We began with the atom in its stationary mixed state (de- 
pending on f2) and the prior distribution of was that 
appropriate to an atom at a random point in a stand- 
ing wave with a maximum Rabi frequency f^max = IO7. 
We analyzed five different measurement schemes, direct 
detection, a particular adaptive scheme [9], homodyne 
detection of the x quadrature, homodyne of the y, and 
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heterodyne. We can summarize the results of the paper 
using four different measures of the effectiveness of the 
measurement. The first two relate to the knowledge ob- 
tained about ri. One is A/;, the long-time (t >> 7~^) 
increase in the average information about the parameter 
il. The other is Vi, the long-time average variance in 
17. The next two relate to the knowledge obtained about 
the system. One is pi, the long-time purity. This mea- 
sures how much is known about the system, given the 
long-time knowledge about the unknown parameter Q. 
The other is Ps, the short-time = a few 7^^) purity. 
This time is long enough that, if Q were known, the sys- 
tem would have been more-or-less completely purified, 
but short enough that the actual amount of information 
obtained about is small. That is, it measures how well 
the measurement can purify the state despite the large 
initial uncertainty in the dynamics. 

The results of our work is summarized in the table 
below, using the four measures of effectiveness for the five 
different detection schemes. Rather than quote figures for 
these four measures, we use a rating system (★ to ★★★★), 
the details of which are explained in the caption. This 
allows the results to be taken in at a glance. 





Detection Schemes 


Measure 


Direct 


Adapt 


Homo X 


Homo y 


Hetero 


Ah 




* 




-k-k-k-k 




Vi 


* 






kkkk 




Pi 






★ 


kkkk 




Ps 


* 






k 





TABLE I: Ratings for the five different detection schemes, for 
four different measures. Four ★s is the best rating and one 
★ the worst. For A/;, any rating above ★ indicates that the 
information about Q, continues to increase with time, with 
the lower cut-offs for and *^*r** being AI; = 2.5 and 5 bits 
respectively att — 507^^. For Vi, any rating above * indicates 
a variance in fl that decreases, with the upper cut-offs for ★** 
and ★★★★ being Vi — 7^ and 107^ respectively sd, t — 507^^. 
For pi, a rating above indicates a purity that continues to 
increase with time. For schemes where the purity saturates, 
the lower cut-off for is pi = 0.9. For schemes where the 
purity continues to increase, the lower cut-off for **** is pi = 
0.95 at t — 507^^. Finally, for ps, the lower cut-offs for 
and *^*r** are, respectively, Ps = 0.65, 0.75, 0.85 at 
t = 87" In all cases fimax = IO7. 

From the table it is observed that homodyne y 
(Sec. IV D) was the best detection scheme by all mea- 
sures except for the short-time purification, for which it 
was the worst. Both of these aspects are explained by the 
fact that this scheme measures cry , the dynamics of which 
depend strongly on CI. Hence the measurement record 
contains a lot of information about O, including its sign 
(because rotations over the top of the Bloch sphere are 
different from rotations under the bottom). This also 
enables the purity to approach unity as time increases. 
However, for short times, when little information about 
fl has been obtained, a y measurement is actually very 



poor for purifying the state. That is because the mea- 
surement tends to produce states with well-defined values 
of y, and these are states that are very sensitive to the 
rotation around the a;-axis at rate VI. For a poorly known 
fi, this tends to make the system state more mixed, so 
that the purity grows only as the information about f2 
increases. 

After homodyne y detection, the method that pro- 
vided most information about Q was direct detection 
(Sec. IV A). Under direct detection, the count rate is 
proportional to -I- 1, and (like ay), the dynamics of 
az depend strongly upon O, due to the Rabi rotations 
around the x-axis. However, in terms of ctz, rotations 
around the -l-x-axis from the ground state are indistin- 
guishable from rotations around the — a;-axis. Hence the 
measurement cannot distinguish the sign of fl and there 
is no change in the ensemble averaged variance as time 
increases. As a consequence, the purity saturates at a 
low value. The short time purification is poor also, for a 
similar reason to that for homodyne y detection. 

The adaptive detection is almost complementary in 
its qualities to homodyne y detection. As explained in 
Sec. IV B, it yields almost no information about fi, be- 
cause the rate of detections in steady state is independent 
of f2. In particular, it yields no information about the 
sign of rt, so the variance is constant. As a consequence, 
the purity does not approach unity. Nevertheless, it does 
approach a quite high value, of over 1 — 7/(-\/2 r^max), 
which is 0.93 for flmax = IO7. This is because the condi- 
tioned states are, for large f2, asymptotically independent 
of 51, as they approach ax eigenstates. This explains why 
the adaptive scheme gives the best results for short-time 
purification: the conditioned states are almost unaffected 
by the uncertainty in il. 

Homodyne x detection (Sec. IV C) is in many ways 
similar to the adaptive scheme, and this is readily un- 
derstandable since it would be expected to produce con- 
ditioned states tending towards a^ eigenstates. Like 
adaptive (and direct) detection, the sign of f2 is inde- 
terminable so the variance is constant. Hence the final 
purity does not approach unity. Although its asymptotic 
value is not as high as that for adaptive detection, it is 
higher than that for direct detection. This is as expected, 
since the conditioned states, being imperfectly localized 
towards the a;-eigenstates, are still affected by ft. This 
also explains why the initial purification is not quite as 
good as for adaptive detection, and why information con- 
tinues to be gained (albeit slowly) as time increases. 

The final scheme, heterodyne detection (Sec. IV E), is 
most easily understood by viewing it as an equal mixture 
of homodyne x and homodyne y detection, which is in 
fact a completely rigorous viewpoint. All of the ratings 
for heterodyne detection are intermediate between those 
for the two homodyne schemes. 

In conclusion, we have shown that gaining knowledge 
about an unknown dynamical parameter by monitoring 
the system is a quite different phenomenon from gaining 
knowledge about the system itself. We have also dis- 
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tinguishcd different sorts of knowledge acquisition with 
distinct characteristics: for the unknown parameter, in- 
formation gain (in bits) versus reducing the variance; and 
for the system, short-time purity gain versus long-time 
purity gain. The ability to acquire knowledge in these 
various ways is extremely sensitive to the choice of mon- 
itoring scheme (which does not affect the average evolu- 
tion of the system). For the system we investigated, ex- 
plaining the particulars of this sensitivity depends upon 
a detailed understanding of the conditional dynamics of 
the system. Our discoveries may have important impli- 
cations for the suitability of different quantum feedback- 



control techniques [23, 24] in experimental systems with 
unknown dynamical parameters. Another direction for 
future work could be to investigate the effect of realis- 
tic imperfections in the detection schemes on state and 
parameter estimation in open quantum systems. 
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